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We present an ab initio theory of transport quantities of metallic ferromagnets developed in the 
framework of the fully relativistic tight-binding linear muffin-tin orbital method. The approach is 
based on the Kubo-Stfeda formula for the conductivity tensor, on the coherent potential approx- 
C , imation for random alloys, and on the concept of interatomic electron transport. The developed 

^ ■ formalism is applied to pure 3d transition metals (Fe, Co, Ni) and to random Ni-based ferromagnetic 

alloys (Ni-Fe, Ni-Co, Ni-Mn). High values of the anisotropic magnetoresistance (AMR), found for 
\^ , Ni-rich alloys, are explained by a negligible disorder in the majority spin channel while a change of 

CNj ■ the sign of the anomalous Hall effect (AHE) on alloying is interpreted as a band-filling effect without 

a direct relation to the high AMR. The influence of disorder on the AHE in concentrated alloys is 
investigated as well. 

o 

C/3 ■ PACS numbers: 72.10.Bg, 72.15.Gd, 75.47.Np 

I. INTRODUCTION 

■ Two galvanomagnetic phenomena discovered in the 19th century, namely, the anisotropic magnetoresistance 
(AMR)^ and the anomalous Hall effect (AHE)^ in ferromagnets, attract ongoing interest both in basic and applied 

I ■ physics. The AMR of bulk systems together with the giant magnetoresistance of magnetic multilayers plays an 
' ^ , important role in magnetic data storage;^ whereas the AHE is not fully understood at present despite the tremendous 
g ■ past and recent research activity.^ However, the basic origin of both phenomena is well known for a long time and 
^ ] it was identified with the simultaneous presence of spin polarization and spin-orbit (SO) interaction. Among open 
I— I- problems, one can list, e.g., correlations between the anomalous Hall conductivities and the longitudinal conductivities 
reported for a number of systems?^ or the occurrence of large AMR values and the change of sign of the AHE at 
^ ■ same composition (around 15 at. % Fe) of random NiFe alloysi^ Let us mention that the AMR and AHE are 

■ related, respectively, to the symmetric and antisymmetric parts of the resistivity tensor of the solid. For this reason, 
, an internally consistent theory of all transport phenomena (resistivity, AMR, AHE, etc.) represents an important 

■ task of condensed-matter physics. 

\ Existing theoretical techniques for investigation of transport properties include the Boltzmann equation and the 
I • Kubo linear response theory which represent appropriate tools to describe electron scattering on impurities or phonons. 
,—1 ■ The latter approach uses the Kubo-Greenwood formula^ for longitudinal resistivities and the AMR while the Kubo- 
T-H . Stfeda formula^ is a starting point for the AHE.^ However, the AHE contains not only an extrinsic contribution 

■ due to the electron scattering, but also an intrinsic one that is solely due to the band structure of an ideal crystal. 
*^ , The intrinsic AHE seems to dominate over the extrinsic part in a number of cases; its values are related to the 

Berry curvatures of the Bloch states of the crystali^ These theoretical concepts and techniques have recently been 
used to address various problems on a model level, such as, e.g., to describe quantitatively the AHE and AMR in 
diluted magnetic semiconductors fi2r— or to develop a unified scheme for the AHE of systems with high and low 
conductivities 

Materials-specific ah initio theory of galvanomagnetic phenomena is a straightforward extension of the model- 
level techniques only for the intrinsic AHE of pure magnetic crystals. ^^"^^^ Transport properties due to impurity 
scattering, in particular residual resistivities, of substitutionally disordered alloys without the SO interaction have 
been systematically studied by means of the Korringa-Kohn-Rostoker (KKR) method and the coherent potential 
approximation (CPA).^^ An alternative first-principles technique has recently been developed within the tight-binding 
(TB) linear muffin-tin orbital (LMTO) method.^" Both approaches agree quantitatively for resistivities of metallic 
alloysSS and diluted magnetic semiconductors The AMR of random cubic alloys - based on diagonal elements 
of the conductivity tensor, i.e., on the Kubo-Greenwood formula - have been studied in the fully relativistic KKR 
method)^"— similar results of the TB-LMTO method with the SO interaction included^^ compare well to those of 
the KKR method again. These studies prove that the SO interaction can have a dramatic effect also on residual 
resistivities of ferromagnetic alloys of light (Sd) elements. 
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An early attempt^ to calculate the AHE in random alloys from the Kubo-Greenwood formula using the KKR-CPA 
method was shown to be incorrect a correct formulation based on the Kubo-Stfeda formula has appeared only very 
recently!^ The authors of Ref. [H have suggested to interpret the coherent part of the anomalous Hall conductivity in 
a random alloy as the intrinsic contribution to the AHE while the incoherent part - the co-called vertex corrections - 
has been identified with the extrinsic contribution to the AHE. Moreover, they have studied the case of dilute alloys 
and have shown that the intrinsic contribution exhibits a well-defined limit for vanishing concentration of impurities 
in contrast to a divergence of the extrinsic contribution, in reasonable agreement with experimental data for FePd 
and NiPd random alloys. 

The purpose of the present paper is to formulate the full conductivity tensor of ferromagnetic metals and substitu- 
tionally disordered alloys in the relativistic TB-LMTO-CPA method in the atomic sphere approximation (ASA) and 
to illustrate its applicability to the galvanomagnetic phenomena of systems containing 3d transition-metal elements 
(Fe, Co, Ni, Mn). The paper is organized as follows. The developed method is presented in Section [Til Section fH Al 
summarizes the most important relations of the fully relativistic TB-LMTO formalism relevant for the subsequent 
development of the transport theory, which is given in Section IH Bl Technical parts of the derivations are left to 
Appendices while details of the numerical procedures employed can be found in Section [H CI The calculated results 
for selected systems and their discussion are contained in Section IHII The AMR of Ni-based alloys is discussed in 
Section UlI Al while Sections IHI Bl and IHI CI are devoted to the AHE in pure metals and in random alloys, respectively. 
The main conclusions of the paper are summarized in Section IIVI 



II. METHOD 



A. Relativistic TB-LMTO-ASA method 



The Hamiltonian in the orthogonal LMTO representation for a ferromagnetic system treated in the fully relativistic 
LMTO-ASA method can be written as^"^ 

i/ = C+(VA)+5°(l-7^°)"'^, (1) 

where the C, a/A and 7 are site-diagonal matrices of potential parameters and the iS" denotes the matrix of canonical 
structure constants. The form of the Hamiltonian H is similar to the non-relativistic or non-magnetic cases^Sr— 
however, the structure of H is more complicated in the spin-polarized relativistic case for two reasons.^^— First, 
the matrices involved in Eq. ([T]) have different kinds of indices: the site index R is combined either with the usual 
relativistic index A — (k/i), or with a composed index A — (ifiX). Here the quantum number /i is related to z- 
component of the total angular momentum, the £ is related to the orbital angular momentum, the non-zero integer k 
is related in a well-known manner to the total angular momentum quantum number j (2j + 1 = 2\k\) and to the £ 
{2£ + l = |2k-|-1|), and the index A labels inequivalent regular solutions of the single-site problem in each ^/i-channel.— 
The Hamiltonian matrix has thus indices H^,ji^, ^^^^ — Hfi'tfiixiju^iX while the structure constant matrix has indices 

•^R'A' RA ~ "^R'kV Rk/^- Second, the site-diagonal matrices C, -s/A and 7 are not fully diagonal in the corresponding 
internal indices, (k/x) or {£iiX), see Refs. [27l - l29l for details. Moreover, the \fK is defined as a matrix in the mixed 
indices (k' fi' , £fiX) which cannot be understood as a square root of any matrix A. However, for the transport theory 
developed in Section IH Bl the detailed structure of the matrices is less important and only their general properties 
are relevant, such as, e.g., the hermiticity of C, 7 and S^. 

The treatment of transport properties of disordered systems requires the Green's functions (GF). The basic one - 
called the physical GF - is defined as the resolvent of the Hamiltonian H ([Ij: 

G{z) = {z-H)-\ (2) 

where z denotes a complex energy variable. Other GF's are introduced in the TB-LMTO formalism.— i^i^^ The 
matrix of screened structure constants S*" in the TB-LMTO representation (superscript a) is defined by 

5" = (3) 

where a denotes a site-diagonal matrix of screening constants, and the site-diagonal matrix of screened potential 
functions P"(z) is defined as22i2S 

P"{z)=\VA{z^C)-\VA)+ +-f-aY\ (4) 
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The auxiliary GF in the TB-LMTO representation is then defined as 

5"(z) = [P"(z)-5"]-\ 



(5) 



which represents a simpler quantity for theoretical and numerical treatments than the physical GF G{z), Eq. 
Both GF's are related to each other by linear rescaling 



G(z) = A"(z) + Ai"(z)5"(z)/i"(z), 

where the quantities A"(z), fJ-^iz) and Il°'{z) denote site-diagonal matrices 

-1 



(6) 



M"(z)-(VArMl + (a-7)P"W] 
m"(z) = [1 + P"(z)(«-7)] 



(7) 



The relation ([6]) is indispensable for an efficient treatment of bulk and layered systems;-^ its proof is sketched in 
Appendix O 

The substitutional disorder in a system on a fixed, non-random lattice can be best treated in the CPA applied to 
the auxiliary GF ([5]). The configurationally averaged g°'{z) is then given by2Si2ii2^ 



(8) 



where the V°'{z) denotes the site-diagonal matrix of the coherent potential functions. Their determination from 
the CPA selfconsistency condition in the relativistic spin-polarized case as well as details of the LSDA selfconsistent 
procedure can be found elsewherci^i^ 



B. Full conductivity tensor 

The conductivity tensor at zero temperature is given according to the Kubo-Stfeda formula-'- as 

a,,, - ao Tr {V,, (G+ - G_) KG_ - V^,G+V, (G+ - G_) 

+ i(X^K-X.K')(G+ -G-)}, (9) 

where the subscripts /i and v denote indices of Cartesian coordinates {fJ,,i^ G {x,y,z}), the trace (Tr) is taken over 
all orbitals of the system, the energy argument of the GF's is omitted since it equals the Fermi energy Ep, and we 
abbreviated G± = G{E-p ± iO). The symbols X^^ and denote, respectively, the coordinate and velocity operators. 
The numerical prefactor do reflects the units employed and the size of the sample; with Ti = 1 assumed here, it is 
given by ctq — / {A-kVqN) , where Vb is the volume of the primitive cell and N is the number of cells in a big finite 
crystal with periodic boundary conditions. 

In analogy to the non-relativistic formulation of electron transport;^ the coordinate operator is represented - 
in an orthonormal LMTO basis leading to the Hamiltonian ([T]) - by a matrix diagonal in the (RA)-index as 

(^A')r'A',RA = -^R-'R -^A'A ^R' (10) 

where Xf^ is the ^-th component of the position vector R. The velocity operator VJ^ is then defined as a quantum- 
mechanical time derivative of AT^ {h = 1): 

V^ = -i[X^,H], (11) 

where [A, B] = AB — BA denotes a commutator. The physical idea behind the simple rule (jlOp is an approximation 
of the true continuous coordinate by its step-like "integer" part that is constant inside each atomic sphere. This 
leads to a systematic neglect of any intraatomic currents so that the resulting conductivity a^i/ describes only the 
net electron motion between neighboring atomic sites. The final simple result for the X^ (jlOp is then obtained 
from basic properties of the phi orbitals \<t'-Eitd phi-dot orbitals used in the definition of the orthonormal 

LMTO basisi^"— In particular, one employs their orthonormality, ('/'r'A' |0ra) = <5r'r^a'A ('^R'A' I'^ra) = 0, 
together with a neglect of small quantities {4>ft.K'\'t''R,id which is consistent with the second-order accuracy of the 
LMTO Hamiltonian (HJ. 
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For practical calculations, one can recast the original form of the conductivity tensor cr^^ ([9]) into a more suitable 
version in which the velocities V"^ ([TTt are replaced by effective velocities defined by 

v'^^-i[X^,S"], (12) 

and in which the physical GF's G± are replaced by the auxiliary GF's g'^ = g°'{EY ± iO), see Eq. ([5]). This transfor- 
mation rests on two relations. The first one connects both velocity operators by 

= {VA)+{F'^r^ v"^ [(F")+]-iyA, (13) 

where 



F" = 1 + S"{a~-f) (14) 

and, consequently, (F")+ = 1 + (a — jjS". The velocity relation (fT3| can easily be obtained from the explicit form of 
H ^ and from an identity S'^{l — jS'^)^^ — {F°')~^S°' ~ ^"[(F")^]^^, which is a direct consequence of the screening 
transformation ([3]). Note that the coordinate operator in the definition of the effective velocity v'j^ (|12p has to be 
understood as a diagonal matrix in the (RA)-index, i.e., 

(^m)r'A',ra = Sr'R Sa'A ^r, (15) 

and that both operators X^, Eqs. pH)) and (|15p . commute with all other site-diagonal operators (C, \/A, 7, a). The 
second relation connects the physical and auxiliary GF's by 

G{z) = iVAy\F-)+ [(a - 7) + 9"{z)F-] [{VA)+]-\ (16) 

which represents a complementary relation to Eq. (j6|) and which is proved in a similar way in Appendix El The use 
of Eqs. ([T3|) and (fT6|l in the conductivity tensor ([9|) yields another form of the latter, namely, 

a^. = ao Tr {z;^ (g^ - 5^) <g2 - <g?< (ff? - g^) 

+ i (XX ~^.«^) (5? -<?-)}• (17) 

The proof of equivalence of both expressions for the a^i, is given in Appendix IbI 

The similarity of both conductivity formulas, Eqs. © and p7)) . is obvious. However, the latter one is better suited 
for configuration averaging since the effective velocities (|12p are non-random operators in contrast to the original 
velocities ([TT]) and to the velocities in the KKR-CPA method. The CPA-average of the cr^i/ (fT7|) can thus be reduced 
to averages of the auxiliary GF's ([8|) and to the standard form of vertex corrections for the CPA-average of products 
g°'{zi)(^g°'{z2)^ The same transformation was achieved in our previous study^" for longitudinal conductivities within 
the Kubo-Greenwood formula^ and the non-relativistic TB-LMTO method. The derived result, Eq. pT)) . represents 
an extension in two directions: to spin-polarized relativistic systems and to the full conductivity tensor, including the 
anomalous Hall conductivities. 



C. Implementation and numerical details 

The fully relativistic theory developed in Section III Bl has been applied to random binary alloys of 3d transition 
metals (Mn, Fe, Co, Ni), where the effects of the SO interaction are rather weak as compared to the alloy bandwidths 
and exchange splittings. For this reason, most of the results shown in the next section were calculated by means 
of a simplified version of the theory, in which the SO coupling was included as an on-site perturbation term of the 
• S form to the LMTO Hamiltonian in the scalar-relativistic approximation (SRA))^!^ This SRA-I-SO approach 
proved to yield results in surprisingly good quantitative agreement with results of fully relativistic LMTO and KKR 
techniques, both for magnetic moments and densities of states^? and for longitudinal resistivities^^ of 3d transition- 
metal systems. In a few selected cases, results of the fully relativistic Dirac (FRD) theorj^ will be shown for a direct 
comparison of both approaches. 

The employed LMTO basis comprised s, p, and d-type orbitals, and the LSDA selfconsistency was achieved with 
the Vosko-Wilk-Nusair exchange-correlation potential^^ and for magnetization pointing along z-axis of the alloy cubic 
structures. More details about the ground-state electronic structure calculations were published elsewhere . ^^1^^''^'' The 
numerical prefactor ctq used in the final conductivity formula, Eq. (fTT]). is given by ao = e^/(47r?iVbiV), where N is 
the number of k vectors sampling the whole Brillouin zone (BZ).^*^ The terms bilinear in the auxiliary GF's g± in 
Eq. (jl7p were configurationally averaged in the CPA with the corresponding vertex corrections included^ whereas 
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FIG. 1. (Color online) The anomalous Hall conductivity a^y of fee Ni versus the imaginary part of energy e for different numbers 
A'' of k vectors in the full BZ. The converged value obtained by methods based on the Berry curvature is axy = 22.0 x 10* S/m, 
see Ref.iia and Table |ll] below. 

the term linear in the g± has been omitted here for symmetry reasons^ (it vanishes identically due to the inversion 
symmetry of the Bravais lattice). In a general case, the latter term in Eq. (|17l) is essentially equivalent to the so-called 
Fermi-sea contribution to the AHE, which can be neglected in metallic systems."'*' 

Particular attention was paid to convergence behavior of the conductivities with respect to the number N of 
sampling k points. We found that large TV's were necessary and that energy arguments of the GF's had to contain 
a non-zero imaginary part, z — Ep ± ie. The quantity e should be as small as possible for purely theoretical values 
of cr^i/; however, for a comparison with experiment, finite values of s can also be used in order to simulate additional 
scattering on phonons, magnons and defects in real samples. Figure [T] displays the anomalous Hall conductivity cTxy 
of fee Ni as a function of e and N. One can see that converged values of (Jxy for small values of e require at least 
N ^ 10^ sampling points in the full BZ. In this study, we used e = 10^^ Ry and « 4 x 10^ for bcc Fe while 
N Ri 9 X 10* were used for fee Co and Ni-based systems. 



III. RESULTS AND DISCUSSION 

For cubic systems with magnetization along z-axis, the conductivity tensor has three independent non-zero matrix 
elements, namely, axx = cTyy, a^z, and axy = —cTyx] the same structure is found for the resistivity tensor p^^. The 
AHE is obtained quantitatively from the non-diagonal elements, axy or pxy, while the diagonal elements define the 
isotropic resistivity p — {2pxx + Pzz)/i and the AMR ratio ^p/ p, where Ap = pzz — Pxx- 



A. AMR in random fee Ni-based alloys 

A brief summary of the present TB-LMTO results and their comparison to available experimental data is shown 
in Table |T] for three Ni-based ferromagnetic alloys with the same average number of electrons per atom Z = 27.7. 
One can see that the resistivities of the Ni-Co and Ni-Fe alloys are smaller by a factor of two to four as compared to 
experimental values.— Similar underestimation of the resistivity values has been obtained by the relativistic KKR- 
CPA method^'^-i while the scattering theory combined with a supercell TB-LMTO approach^ yields resistivities that 
are only ^ 20% smaller than the measured ones. The calculated AMR values of Ni-Co and Ni-Fe values exceed slightly 
the measured ones while the opposite discrepancy is found for the ferromagnetic Ni-Mn system: for a Nio.97Mno.03 
alloy, the measured AMR amounts to 7.3%^^ whereas the present calculation yields the AMR of 3.3% and 2.4% in 
the FRD and the SRA-l-SO approach, respectively. 

The first application of the developed scheme has been done for the AMR of fee Ni-based alloys. A preliminary 
study of the Ni-Co and Ni-Cu alloy systems including a comparison with the results of the KKR method was published 
elsewhere?^ here we show the AMR for the fee Nii.j-Fe^; alloy calculated in the SRA-l-SO approach, see Fig. [21 One 
can see relatively high values of the AMR for the Ni-rich alloys with the maximum AMR ratio slightly exceeding 20%. 
Similar trends and values were obtained in the fully relativistic KKR method§.'.iS and in experiments;^ 
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TABLE I. Comparison of the calculated and measured transport quantities for three fee Ni-based alloys with the average number 
of electrons per atom Z = 27.7: the isotropic resistivity p (in /^f^cm), the AMR ratio and the anomalous Hall conductivity G^y 
(in kS/m). The displayed calculated values correspond to the FRD approach; in parentheses, results of the SRA+SO approach 
are shown. The experimental values are taken from Ref. [i^ . 



alloy 


P 


AMR 




NIo.tCoo.s calc. 


1.06(0.70) 


0.47(0.42) 


-88(-74) 


exp. 


2.8 


0.3 




Nio.85Feo.15 calc. 


2.21(1.57) 


0.29(0.24) 


-40(-29) 


exp. 


4.2 


0.18 




Nio.gMuo.i calc. 


22.8(23.6) 


0.009(0.004) 


65(54) 



0.3 



0.2 



Q. 



0.1 




0.3 

Fe concentration 



0.6 



FIG. 2. (Color online) The calculated values of the anisotropic magnetoresistance of random fee Nii_a;Fea; alloy as a function 
of Fe concentration and their comparison to experimental values.— 



The high AMR values for the Ni-rich Ni-Fe alloys and even higher theoretical^^^ as well as experimental^!^ values 
for the Ni-rich Ni-Co alloys deserve attention and call for explanation. Existing theoretical approaches to the AMR 
in metals and alloys are based on inclusion of the SO interaction in the lowest-order perturbation expansion42iiliM 
However, the validity of such schemes for the Ni-rich alloys is limited due to the well-known fact that the SO interaction 
increases the resistivity by a factor of two or more as compared to the resistivity within the two current mode h^^i^^ 
In order to reveal the possible origin of the high AMR in these systems, we considered also the fee Ni-Mn alloy 
with a ferromagnetic order in the whole concentration range studied (which represents the true ground state only for 
Mn concentrations below 15 at.% M n^^'^° ) and an artificial fee Ni-Fe(*) alloy in which the majority spin (spin-up) 
potential of Fe atoms was replaced by that of Ni atoms leading thus to a system with no spin-up disorder. The real 
Ni-Co and Ni-Fe alloys are featured by a very weak disorder in the majority spin channel)^ the four ferromagnetic 
alloys considered in the present study form a set in which the spin-up disorder decreases along the sequence Ni-Mn, 
Ni-Fe, Ni-Co, and Ni-Fe(*). 

The resulting AMR values of the four systems are plotted in Fig. [3] as functions of the change of the number of 
valence electrons AZ due to alloying. The latter quantity is defined as AZ(a;) — x for the Nii_a;Coa; alloy, AZ(a:) = 2x 
for the Nii_a;Fej; alloy, and AZ{x) — 3x for the Nii-a^Mn^; alloy. The calculated AMR increases along the sequence 
Ni-Mn, Ni-Fe, Ni-Co, Ni-Fe(*), in qualitative agreement with experiment for the real alloys4^ The obtained trend 
indicates that the AMR ratios are directly related to the disorder strength in the majority spin channel, with the 
highest AMR ratios obtained for alloys with negligible spin-up scattering. This conclusion does not depend on the 
particular approach (SRA+SO, FRD) employed in the calculation, as can be seen from a comparison of results for 
AZ = 0.3, see Table III 
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FIG. 3. (Color online) The calculated values of the anisotropic magnetoresistance of random fee Ni-based alloys as functions 
of the valence charge difference AZ. 



TABLE 11. The calculated and experimental values of the anomalous Hall conductivity Gxy (in kS/m) for 3d transition-metal 
ferromagnets. The displayed values obtained in this work correspond to the FRD approach; in parentheses, results of the 
SRA+SO approach are shown. 





bcc Fe 


fee Co 


fee Ni 


this work 


-65(-55) 


-36(-34) 


241(243) 


Berry curvature 


-75" 


-25'' 


220*= 


KKR method 


-64^* 




164'* 


experiment 


-103*= 







" Reference llQ. Reference ll§. Reference |17|. ** Reference i26i. *= Reference 1 5 ll. ^ Reference (52^ 



B. AHE in elemental ferromagnets 

The calculated anomalous Hall conductivities axy of three cubic 3(i transition metals are listed in Table [IT] together 
with other theoretical and experimental results. The Co was treated here in the fee structure, which represents the 
natural cubic approximation to the ground-state hexagonal close-packed structure.^'* It should be noted that our sign 
convention for the axy differs from that adopted in the previous studies all signs in the table are thus taken 
consistently with our convention. 

One can see that the results of the TB-LMTO method agree reasonably well with the data based on the Berry 
curvatur&iSrlS, and on the fully relativistic KKR Green's-function technique. The biggest discrepancy in Table HIl 
related to the measured and calculated AHE of fee Ni, might reflect effects of strong correlations on the electronic 
structure 1^ 



C. AHE in random fee Ni-based alloys 

The effect of alloying on the AHE has been studied for the random fee Ni-based alloys discussed in Section IIII Al 
In order to assess the role of chemical disorder on the AHE, the W\i-xCox alloy has also been treated in a simple 
virtual crystal approximation (VGA) in which the alloy constituents are replaced by an effective atom with the atomic 
number Z{x) = 2^ — x\ the resulting system is denoted as Ni-Co(*). 

The concentration trends of the axy, calculated in the SRA-I-SO approach, are summarized in Fig. We note that 
the dilute limit of the alloys, where the axy diverges due to the diverging incoherent contribution,^^ has not been 
studied here. All studied systems exhibit a change of the sign of the AHE due to alloying. This can be understood as 
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FIG. 4. (Color online) The calculated values of the anomalous Hall conductivity a^y of random fee Ni-based alloys as functions 
of the valence charge difference AZ. 



a consequence of band filling and of the opposite signs of the a^y for pure fee Co and Ni, see Table |IT1 A similar sign 
change has been observed a long time ago for Ni-rich Ni-Fe and Ni-Co alloys i^ i^^'^^ This feature was related to the 
occurrence of the high AMR values encountered in these systems with roughly the same composition.— However, 
our calculated data for the ferromagnetic Ni-Mn alloy witness, that the sign change of the AHE can be obtained also 
for a system with much smaller AMR values. 

Let us discuss finally the role of disorder on the AHE from a comparison of the calculated axy for the Ni-Fe and 
Ni-Fe(*) systems as well as for the Ni-Co and Ni-Co(*) alloys. As one can see from Fig. SI the neglect of the disorder 
in the majority spin channel (Ni-Fe) or in both spin channels (Ni-Co) has a strong influence on the resulting a^y- 
This analysis indicates that results of simplified treatments of alloying, such as, e.g., in Ref. [s^, should be taken with 
caution. Moreover, the o-^y^ for the Ni-Co(*) alloy differs significantly both from the total a^y for the corresponding 
Ni-Co alloy in the CPA and from its coherent part, cr™'^. For the equiconcentration Ni-Co alloy, we get: cr^y^ — ~14 
kS/m, (Txy = —88 kS/m, and u™'' — —93 kS/m. This example proves that the coherent part of the anomalous 
Hall conductivity of a real concentrated alloy is not related directly to that of the corresponding effective crystal, in 
contrast to the case of the dilute alloys 



We have reformulated the Kubo-Stfeda expression for the conductivity tensor—^ in the fully relativistic TB-LMTO- 
CPA theoryf22. This extends the applicability of our previous transport formalism^^ to all elements of the conductivity 
tensor, indispensable for the anomalous Hall effect. The effect of the spin-orbit interaction has been included and 
implemented numerically in the Dirac four-component formalism as well as in a simple perturbative manner. The for- 
mer approach is inevitable for systems with heavy elements while the latter scheme is advantageous for interpretation 
of the results in terms of the underlying electronic structure owing to the non-relativistic (or scalar-relativistic) basis 
set of orbitals. First applications to galvanomagnetie phenomena in pure metals (Fe, Co, Ni) and in ferromagnetic 
Ni-based alloys yield results in reasonable agreement with other techniques. The performed calculations for the real 
and artificial alloy systems clarified some aspects of the anisotropic magnetoresistance and the anomalous Hall effect 
and their possible interrelation. The developed technique has already been generalized and applied to alloy structures 
with several sublattices and different degree of atomic ordering. Further questions, such as, e.g., the influence of com- 
plex magnetic orders on the transport properties in systems with competing ferro- and antiferromagnetic interactions, 
remain a task for future. 



IV. CONCLUSIONS 
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Appendix A: Relations between the Green's functions 



The proof of Eqs. (jG]) and (|16p is based on a general identity valid for matrices x, x, y, y, and /, that are coupled 

by 

x^x{l + fx)-^ = (1 + xfy^x, x^x{l- fx)-^ = (1 - xfy^x, 

y^yii + fyr' = ii + yfr'y, y ^y{i- fy'r' = i^-yfr'y, (Ai) 

if corresponding inverses exist. This yields relations 

x-y={l + xf)-\x-y){l + fy)-\ 
(x - y)-' = (1 + fy){x - + xf), (A2) 

and the general identity 

(x-y)-' = -{l + fx)f + {l + fx){x-y)-\l+xf). (A3) 

This identity can be used for x = [{y/A)+]-\z - C){VA)-\ y = S'°(l-75°)"\ and / = q;-7, so that x = P"(z) 
and y — S". Substitution into Eq. (jA3l) and subsequent multiplication by (a/A)^^ from the left and by [(-\/A)~'"]~^ 
from the right leads to the relation ([5]) with the matrices A"(z), i^"{z) and p."{z) given by Eq. ([7]). Similarly, the 
identity dMl) can be applied to y = [{VA)+]-\z - C){^/A)-\ x = S°{l-jS°)-^, and f ^a~-f, so that y = P"(z) 
and X — S". This yields the complementary relation between the Green's functions, Eq. ([2 



Appendix B: Equivalence of expressions for the conductivity tensor 

For the proof of equivalence of Eqs. 1^ and p7)) . we drop the superscript a at v'^, g'^, and F", and abbreviate 
a = a — 7. The definition (fT4|) is thus written as 



F = l + Sa, (Bl) 
and the two relations (fT3|) and ()16p are now written as 

= (\/A)+p-it;^(F+)-iVA, G± -(v^)-ip+(a + 5±P)[(VA)+]-i. (B2) 
Substitution of Eq. (jB2p into the original Eq. © yields after trivial modifications 

crf,„ = (Jo Tr {P^^u^ (5+ - 9-)v^ {a + g^F) - (a + g+F) F~'^v^ {g+ - g^) 

+ i {X,,F-'v, - X.p-ii;^) (5+ - F} = + aj^ , (B3) 

where the a^^J and ct^^? comprise, respectively, all terms bilinear and linear in g±. The first contribution is obvious, 

'^i^J = f^o Tr {Vf, {g+ - .g_) v^,g^ - v^,g+v^ {g+ - ff-)} , (B4) 

while the second contribution is more complicated. It has a form 

r7(2j =aoTr{7V^„(.g+-.g_)}, (B5) 

where the N^^ can be written in terms of the effective velocities (IT^ as 

N^, = v,aF-^v^ - v^aF-\, + iF {X^F-\, - X^F'^v^) 

= i ([X^, 5]a + FX^) P-^I;, - i ([X„ 5]a + FX,) F-\^. (B6) 

Here, the last two brackets can be modified using Eq. (|Bip and [A"p,a] — 0, so that [Xp,F] — [X^,S]a, hence 
[X^, S'Ja + FXf, = X^F and Eq. (IB6| reduces to 

7V^, = i(X^«,-X,w^). (B7) 
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Consequently, the second contribution in Eq. (|B3I) simplifies to 

of, = ao Tr {i {X^v, ~ X,v^) {g+ - g^)} . (B8) 
The sum of Eq. ()B4p and Eq. (|B8|) is now identical to the transformed expression for the a^i, ([T7|) . 
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